Towards an automatic ordering for a symmetric sparse direct solver

نویسندگان

  • Iain S. Duff
  • Jennifer A. Scott
چکیده

In recent years, nested dissection has grown in popularity as the method of choice for computing a pivot sequence for use with a sparse direct symmetric solver. This is particularly true for very large problems. For smaller problems, minimum degree based algorithms often produce orderings that lead to sparser matrix factors. Furthermore, minimum degree orderings are frequently significantly cheaper to compute than nested dissection. In this report, we look at whether we can predict which ordering will be better, using only the sparsity pattern of the matrix. Our aim is to choose efficiently a good ordering for a wide range of large problems from different application areas.

برای دانلود متن کامل این مقاله و بیش از 32 میلیون مقاله دیگر ابتدا ثبت نام کنید

ثبت نام

اگر عضو سایت هستید لطفا وارد حساب کاربری خود شوید

منابع مشابه

Effective Preconditioning through Ordering Interleaved with Incomplete Factorization

Consider the solution of a sparse linear system Ax = b when the matrix A is symmetric and positive definite. A typical iterative solver is obtained by using the method of Conjugate Gradients (CG) [15] preconditioned with an incomplete Cholesky (IC) factor L̂ [4]. The latter is an approximation to the (complete) Cholesky factor L, where A = LL . Consequently, the process of computing L̂ relies to ...

متن کامل

Performance of a Fully Parallel Sparse Solver

The performance of a fully parallel direct solver for large sparse symmetric positive deenite systems of linear equations is demonstrated. The solver is designed for distributed-memory, message-passing parallel computer systems. All phases of the computation, including symbolic processing as well as numeric factorization and triangular solution, are performed in parallel. A parallel Cartesian n...

متن کامل

A Large-Grain Parallel Sparse System Solver

The eeciency of solving sparse linear systems on parallel processors and more complex multicluster architectures such as Cedar is greatly enhanced if relatively large grain computational tasks can be assigned to each cluster or processor. The ordering of a system into a bordered block upper triangular form facilitates a reasonable large-grain partitioning. A new algorithm which produces this fo...

متن کامل

The design and use of a sparse direct solver for skew symmetric matrices

We consider the direct solution of sparse skew symmetric matrices. We see that the pivoting strategies are similar, but simpler, to those used in the factorization of sparse symmetric indefinite matrices, and we briefly describe the algorithms used in a forthcoming direct code based on multifrontal techniques for the factorization of real skew symmetric matrices. We show how this factorization ...

متن کامل

PSPASES: An Efficient and Scalable Parallel Sparse Direct Solver

Many problems in engineering and scienti c domains require solving large sparse systems of linear equations, as a computationally intensive step towards the nal solution. It has long been a challenge to develop e cient parallel formulations of sparse direct solvers due to several di erent complex steps involved in the process. In this paper, we describe PSPASES, one of the rst e cient, portable...

متن کامل

ذخیره در منابع من


  با ذخیره ی این منبع در منابع من، دسترسی به آن را برای استفاده های بعدی آسان تر کنید

برای دانلود متن کامل این مقاله و بیش از 32 میلیون مقاله دیگر ابتدا ثبت نام کنید

ثبت نام

اگر عضو سایت هستید لطفا وارد حساب کاربری خود شوید

عنوان ژورنال:

دوره   شماره 

صفحات  -

تاریخ انتشار 2005